(N 



O 



The instantaneous shear modulus in the shoving model 

Jeppe C. DyrCl 

DNRF Centre "Glass and Time", IMFUFA, Department of Sciences, 
Roskilde University, Postbox 260, DK-4000 Roskilde, Denmark 

Wei Hua Wan^ 

Institute of Physics, Chinese Academy of Sciences, Beijing 100190, P.O. Box 603, China 

(Dated: March 1, 2013) 

We point out that the instantaneous shear modulus Goo ol the shoving model for the non- Arrhenius 

temperature dependence of viscous liquids' relaxation time is the experimentally accessible high- 

^— V ' frequency plateau modulus, not the ideahzed instantaneous affine shear modulus that cannot be 

fvj , measured. Data for a large selection of metallic glasses are compared to three different versions 

of the shoving model. The original shear-modulus based version shows a slight correlation to the 

j^^' Poisson ratio, which is eliminated by the energy-landscape formulation of the model in which the 

I ^ \ bulk modulus plays a minor role. 

'^ ' I. INTRODUCTION 

In a recent ConnnunicatiorilPuosi and Leporini showed from computer simulations that the relevant high-frequency 

^+H ' shear modulus controlling the relaxation is not the idealized shear modulus corresponding to afhne deformations at 

truly infinite frequency. Rather, it is the shear modulus referring to time scales that are on the one hand much shorter 

than any relaxation time, but on the other hand much longer than typical vibration times. This confirms findings by 

C^ other groupsi^""— Puosi and Leporini further proposed an extension of the shoving model to allovir for heterogeneities 

U ■ and shovifcd that the new model fits simulation data very well. This extension is consistent with previous works by 

I ' Khronik et al., who introduced the idea of a distribution of local shear moduli to explain sub Tg relaxations within 

' O . the shoving-model frameworki^ 

^ ' Commenting on the instantaneous shear modulus of the shoving model Puosi and Leporini wrote; "It is quite ap- 

parent that Goo, the central quantity of the standard elastic models, poorly correlates with the structural relaxation 
time." This conclusion derives from the understanding that the shoving model and related elastic models are based 
\ on the idealized affine infinite-frequency shear modulus, not the experimentally measurable high-frequency plateau 
\^ • modulus, which was traditionally used for comparing model to experiment. Unfortunately, theorists and experimen- 
I ' talists have used for many years the same symbol. Goo, for different physical properties. In experiment the term 
—^ 1 "instantaneous shear modulus" always meant the shear modulus at the highest obtainable frequencies (typically MHz 
\^ or GHz, or even high kHz, depending on the technique used)^"— where it generally becomes frequency independent 
\^ , and is usually - though not always - denoted by Goo '^'^ In his excellent text book The dynamic properties of super- 
_^ cooled liquids from 1976 Harrison refers to this quantity as "the high-frequency limiting shear modulus" .* Likewise, 

("^ ] whenever the Maxwell relaxation time tm — tj/Goo in the literature has been calculated for a liquid. Goo was always 
^sj ' identified as the experimental limiting high-frequency shear modulus of the plateau, much below phonon frequencies. 
In liquid-state theory. Goo was traditionally the idealized, truly infinite- frequency limit of the fiuctuation-dissipation 
theorem expressiouji^ii^ proportional to the mean-square shear-stress equilibrium fluctuation. Previously, there was 
• i-H ' no reason to believe that these two quantities should differ in any significant way. It now turns out that for some 
/\ , systems, as temperature changes at constant density, one quantity increases and the other decreases i^ii^ The purpose 
JH ■ of the present paper is to show that the shoving model derivation assumes Goo is the quantity, which Puosi and 
- - - Leporini referred to as the plateau modulus, ^^"^ We suggest a consistent notation for the two instantaneous shear 
moduli, which can be used whenever there is a risk of confusing themiiiii Finally, we reanalyze data for a large 
selection of metallic glasses in order to compare to predictions of three versions of the shoving model. 

Section II reviews the shoving model. Sec. HI recalls the model energy-barrier calculation and relates it to the 
reversible-work theorem of statistical mechanics. Section IV compares model predictions to data, and Sec. V gives a 
brief summary. 

II. THE SHOVING MODEL 

The shoving model for the non- Arrhenius temperature dependence of the main (alpha) relaxation time r of a glass- 
forming liquid predicts that t{T) = roexp [Goo{T)Vc/{kBT)\M^ Here tq is a prefactor of order 0.1 picosecond. Goo 
is the "instantaneous" shear modulus, Vc is a characteristic volume of order a molecular volume, fcs is Boltzmann's 



o 



constant, and T is the temperature. In order to minimize the number of free parameters, the shoving model assumes 
ad hoc that Vc is temperature independent. In this way the model connects directly two experimentally measurable 
quantities, t{T) and Goo{T). The model does not consider what causes the unusually large temperature dependence 
of glass-forming liquids' instantaneous shear moduli. This problem was addressed recently in interesting papers by 
Brito and Wyart, who proposed that the increase of the instantaneous shear modulus derives from a stiffening of the 
Boson peak as temperature is lowered i^"'^^ 

The shoving model and related elastic models have been confirmed for a number of organic, oxide, chalcogenide, 
oxynitride, and metallic liquids fi^i^^—, as well as in some computer simulations j^^'' but failures of the model have also 
been report ed j^^i'^° Most experimental conformations relate to a liquid's temperature-dependent equilibrium relaxation 
time, but there are also tests confirming the shoving and related elastic models for aging experiments i^i^i^^ More 
data are certainly needed before it is clear whether the shoving and related models account generally for the non- 
Arrhenius temperature dependence of r observed in supercooled liquids. For reviews of elastic models the reader is 
referred to Refs. |3j andlSJl 

The shoving model is based on the assumption that the main contribution to the activation energy for a "flow event" 
- a rearrangement from one potential energy minimum to another - is the work done in shoving aside the surroundings 
in order to increase the volume available for rearranging the molecules. The model assumptions areJ^'^'^ 

• The main contribution to the activation free energy is elastic energy. 

• This elastic energy is located in the surroundings of the rearranging molecules. 

• The elastic energy is shear elastic energy, i.e., not associated with density changes. 

To make things simple, the shoving model assumes spherical symmetry of a flow event and calculates the activation 
free energy as the work done in shoving aside the surroundings (expanding a sphere) in order to create room for a 
flow event. 

Is it reasonable to assume that the main contribution to the activation energy comes from the surroundings? What 
about contributions from the rearranging molecules themselves? A simple argument shows that the former contribution 
dominates^i Suppose that rearranging at constant volume is energetically very costly because the molecules are forced 
into close contact during the rearrangement process; this is the main physical idea of the shoving model. In this case, 
allowing for just a slightly larger volume for the rearranging molecules implies a considerable lowering of the energy 
cost. If the radius change is Ar and the energy barrier contribution from the rearranging molecules within the 
sphere is /(Ar), the total energy barrier involves a further quadratic contribution from deforming elastically the 
surroundings: /S.E = f{Ar) + A{Ar)'^ . The fact that the function /(Ar) decreases signiflcantly when Ar increases 
slightly above zero is expressed mathematically as |dln//dlnAr| ^ 1. Optimizing Ar in order to find the lowest 
barrier leads to /'(Ar) -I- 2AAr — 0. Thus the ratio between the "shoving" work and the "inner" contribution is 
A(Ar)V/(Ar) = -/'(Ar)/[2Ar/(Ar)] = |ciln//dln Ar|/2 > 1. 

It is the high-frequency shear elastic constant that enters into the shoving model prediction because the expansion 
of a sphere in an elastic solid results in a radial displacement in the surroundings, which varies with distance r to 
the sphere center as r~^^ This is a pure shear deformation, i.e., with zero divergence and thus no density changes 
anywhere (compare to the Coulomb electric field of a point charge oc r~^, which also has zero divergence). Less 
idealized geometries would result in some density change and thus also involve the bulk modulus, but it has been 
shown generally that the bulk elastic energy of the far field of an elastic dipole constitutes a most 10% of the total 
elastic energy.'^^ This confirms the interesting phenomenon of "shear dominance" noted some time ago in various 
contexts of condensed matter physics and materials science (see Refs. |3a and 1371 and their references). 

III. CALCULATING THE FREE ENERGY BARRIER USING THE REVERSIBLE- WORK THEOREM 

OF STATISTICAL MECHANICS 

The idealized instantaneous affine shear modulus refers to the hypothetical situation where one imposes an instan- 
taneous, perfectly affine shear deformation on the system, corresponding to a time scale so fast that the atoms do 
not move, i.e., on a femtosecond time scale. This quantity is not experimentally accessible for the following reasons. 
It is difficult to imagine imposing an affine shear deformation on a system on a femtosecond time scale. Even if 
this were possible, one would not observe a high-frequency limiting shear modulus above THz frequencies - inertial 
effects would set in and cause the modulus to go to zero as frequency diverges without limit. On the other hand, in 
a computer simulation an instantaneous affine shear deformation is easily imposed on a system. 

The energy barrier is calculated as the work done in shoving aside the surroundings. According to a fundamental 
theorem of statistical mechanics, the difference in free energy between two states can be calculated as the reversible 



isothermal work done to bring the system from one to the other state. Thus, in calculating the energy barrier, 
the expansion of the sphere must take place so slowly that there is equilibrium in the surroundings throughout the 
expansion process. On the other hand, since the activation (free) energy refers to the thermally activated creation of 
extra volume in the fixed, glassy structure of the surrounding molecules, the expansion must be fast enough that no 
relaxations take place in the surrounding liquid. This means that the time of expansion must be much smaller than 
the alpha relaxation time. In conclusion, the relevant elastic constants of elastic models refer to high frequencies, still 
much below phonon frequencies. In this frequency range the bulk and shear elastic constants are usually frequency 
independent for highly viscous liquids (ignoring possible secondary relaxations). 

Over the years we have occasionally met the misunderstanding that the "shoving" process is the work done during 
the actual barrier transition of the rearranging molecules, which takes place over a few picoseconds. This is not correct; 
the (free) energy barrier is a difference in free energy between two states - the starting state and the transition state - 
and as detailed above this quantity is calculated by reference to statistical mechanics and the reversible- work theorem. 

Both the experimentally measurable high-frequency plateau shear modulus - referring to frequencies much below 
phonon frequencies - and the experimentally non-accessible truly instantaneous afhne shear modulus have traditionally 
been denoted by Goo- This would not present a serious problem if the two moduli were more or less identical or, 
from the shoving model perspective, if they were proportional in their temperature variation throughout the phase 
diagram. It now appears that neither is the case.^^^ This calls for introducing an unambiguous notation. We suggest 
denoting the experimental "instantaneous" shear modulus by Goo,p (p for plateau) and the idealized, affine truly 
instantaneous shear modulus by Goo,af ("af" for affine). In this notation, which need only to be used whenever there 
is a risk of confusion, the shoving model prediction is 



r(T) = roexp 



GooAT)Vc 



knT 



(1) 



Before proceeding we note the fluctuation-dissipation (FD) theorem expressions for the two instantaneous shear 
moduli. Recall that if one defines S^y = J2i ^iFy.i where Xi is the x component of the position vector of particle i and 
Fy^i is the y component on the force on this particle, the FD theorem expression for the frequency-dependent shear 
viscosity rj{u) is (where V is the system volume) 

, , J,^{S,y{0)S^y{t))e-^-'dt 

v{^) = — r^FT^ ■ 2) 

Since the frequency-dependent shear modulus is related to the viscosity by G(w) = iujri{uj), letting frequency go to 
infinity in Eq. ^ gives the well-known expression for the idealized affine infinite- frequency shear modulus Goo.af r^ 

ri — ^ xyl /„N 

C.oo,af - j^^^y {6} 

As mentioned, this quantity cannot be measured in experiment. The quantity that can be measured, the plateau 
modulus Goo,p, is given by the analogous expression where Sxy is averaged over a few molecular vibration periods 
(i.e., over some picoseconds): 

2 

^°°'P - k^fv ■ ^^' 

We finally note that Goo,af/ {pT) and Goo,p/{pT) are both isomorph invariants. The invariance of the first expression 
was shown in Rcf. 38, that of the latter follows by analogous arguments. This shows that the shoving model survives 
the "isomorph filter" according to which any generally applicable model for the non-Arrhcnius temperature dependence 
of a supercooled liquid's relaxation time must give this as a function of an isomorph invarianti^ 

IV. COMPARING THE ENERGY-LANDSCAPE VERSION OF THE SHOVING MODEL TO 

METALLIC GLASS DATA 

The energy-landscape justification of the shoving modeP^ is based on a classical argument, which estimates the 
barrier height for a jump between two (free) energy minima from the curvature at the minima, leading to In r ex 1/ {v?) 
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FIG. 1: Shear and bulk moduli for a range of metallic glasses (open and full symbols, respectively, left axis) and the glass 
transition temperatures Tg (half- full symbols, right axis). These data are used below for testing the different versions of the 
shoving model. 



where (u^) is the vibrational mean-square displacement 42r— If the high-frequency shear and bulk plateau moduli, 
Goo,p and Koc,p, differ from the ideal, afhne moduli Goo^ai and Koc,!i[, the former are the relevant ones for the phonon 
spectra. In fact, one way to probe plateau moduli is to measure the linear (i.e., low-wavevector) parts of the phonon 
spectra or, equivalently, the shear and longitudinal high-frequency sound velocities. This is also how Ribero et al. 
probed Goo in computer simulations.'^''* In this section term "instantaneous" moduli implies the plateau moduli. 

In a simple, isotropic elastic model the vibrational mean-square displacement may be estimated from the instan- 
taneous moduli by using the fact that for a given wavevector, there are two transverse and one longitudinal phonon. 
li M = K + (4/3)G is the longitudinal modulus, this leads-- to {u'^)/T oc l/Moo,p + 2/Goo,p. For the temperature 
dependence of the relaxation time this implies via Inr oc l/(u^) 



t{T) 



Toexp 



V' 



ksT \MooAT) GooAT) 



(5) 



where V^ is a molecular-sized volume. As shown in Ref. |39J, this expression implies shear dominance for the tem- 
perature dependence of t. This was shown by first defining for any quantity Q that increases as T decreases its 
temperature index as Iq = —d\uQ/dlnT. Equation ([5]) implies that the temperature index of r's activation energy 
can be written (1 — a)lG^ + chIk^ , where < a < 0.08 is obeyed no matter what is the ratio Goo,p/Koo,p-^- In 
other words, at least 92% of the non-Arrhenius temperature dependence of the relaxation time derives from Goo.p's 
temperature dependence. The instantaneous bulk modulus i^oo.p plays only a minor role for three reasons: 

• Two phonons arc transversal for each one that is longitudinal. 

• The instantaneous shear modulus Goo.p affects also the longitudinal phonons. 

• Longitudinal phonons are harder than transverse and therefore contribute less than one third to the vibrational 
mean-square displacement. 

It follows from Eq. ^ that the temperature-dependent activation (free) energy AE{T) is given by 



AE{T) K Goo,p(r)K' 



j^oo,p(r) + 4Goo,p(r)/3 

2i^oo,p(r) + llGoo,p(T)/3 ' 



(6) 



Traditionally the shoving model is compared to experiment by plotting the logarithm of the relaxation time (or, 
equivalently, the viscosity) as a function of Goo,p(r)/T''^^^'^— '^ Such a plot checks whether the predicted linear 
relationship between Inr and Goc,p{T)/T is observed with a physically reasonable prefactor. A convenient way 
of making this plot is via a generalized Angell plot in which the x coordinate is normalized to unity at the glass 



transition by defining X = [Goo,p{T)Tg]/[Goo,p{Tg)T]. Unfortunately, measuring the high-frequency plateau shear 
modulus of the equilibrium metastable supercooled liquid is difficult and not too many data are available even today. 
An alternative way of testing the shoving model is by making use of the facts that 1) for a fixed cooling rate the 
relaxation time has a certain value at Tg, 2) the (DC) shear modulus of the glass, G, is almost temperature independent 
and equal to the instantaneous shear modulus of the liquid at Tg: G = Goo,-p{Tg). The glass shear modulus is easy 
to measure, and the prediction of the model is that GVc/ksTg = Const., in which it is reasonable to assume that 
the constant is more or the same for similar systems. For the range of metallic glasses compared below to model 
predictions we further assume that the characteristic volume Vc is proportional to the molar volume Vm with a 
universal proportionality constant. The point of these simplifying assumptions is to eliminate all non-trivial free 
parameters. 

Figure 1 gives bulk and shear modulus data for a range of metallic glasses, as well as their glass transition temper- 
atures ranging from 317K for some of the Ca-based glasses to 930K for some Fe-based ones. Most of these data were 
discussed and compared to elastic model predictions in Refs. I26ll28l but the cupper-based glasses were replaced. 

Figure 2 compares data to different versions of the shoving model [Eq. ([T])]. Since the glass transition takes place 
when the relaxation time upon cooling reaches a certain value, the standard shoving model prediction is 

Tg^AE{Tg)^GVm (7) 

with a universal proportionality constant. Figure 2(a) tests this prediction based on the data of Fig. 1, where the 
metallic glasses are sorted according to their Poisson ratio {R is the gas constant). Given the diversity of the glasses, 
the model's simplicity, and the simplifying assumption that Vc oc Vm with a universal proportionality constant, the 
data show good agreement with the shoving model. There is, as noted in Refs. (261 and |28| . some correlation with 
Poisson's ratio4i Before addressing whether the energy-landscape version of the shoving model rectifies this, we show 
in Fig. 2(b) how data compares to the elastic model in which G is replaced by K. Here the correlation with Poisson's 
ratio is much stronger. Altogether, Figs. 2(a) and 2(b) show that if an elastic model is to work for metallic glasses 
with a universal value of the characteristic volume Vc in terms of the molar volume Vm, the bulk modulus can play 
only a minor role. This constitutes an experimental demonstration of "shear dominance" . 
For the energy-landscape version of the shoving model Eq. ^ translates into the prediction 

T,ocAi.(T,)ocGV.^f±i|^. (8) 

This is compared to experiment in Fig. 2(c). The fit is good and there is no correlation to the Poisson ratio. 

The three "zero-parameter" elastic models of Fig. 2 have standard deviations, respectively, of 8% for the standard 
shoving model [Fig. 2(a)], 22% for the bulk modulus elastic model [Fig. 2(b)], and 8% for the energy-landscape 
version of the shoving model [Fig. 2(c)]. The latter gives no better fit to the data than the original shoving model 
(the standard deviation of it is 0.6% lower, but this is hardly significant), but has the advantage of eliminating the 
correlation to the Poisson ratio. 

A pragmatic one-parameter version of the elastic models is to allow for an arbitrary combination of G and K 
by assuming that Tg is controlled by aG -I- (1 — a)K. Due to shear dominance one expects a to be close to unity. 
Indeed, Ref. [2^ showed that for a — 10/11 a very good fit to data is obtained for metallic glasses, a fit which like 
the energy-landscape shoving model eliminates the correlation to the Poisson ratio. For this one-parameter model the 
standard deviation is 6%, i.e., somewhat better than the standard shoving model and its energy-landscape version. 

V. SUMMARY 

The instantaneous shear modulus of the shoving model refers to the plateau modulus that is traditionally in 
experiment denoted by Goo. This quantity is measurable, in contrast to the idealized affine instantaneous shear 
modulus of liquid state theory, which was also traditionally denoted by Goo ■ A consistent notation has been suggested 
for distinguishing between these two quantities whenever there is a risk of confusing them. Data for a range of metallic 
glasses have been shown to be consistent with shoving model predictions; in particular the energy-landscape version 
of the model eliminates the correlation of model predictions to Poisson's ratio. 
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FIG. 2: Comparing shoving model predictions to data for a range of metallic glasses (Fig. 1). In order to reduce the number of 
free parameters as much as possible it is assumed that the characteristic volumes Vc and Vj via Eqs. ([1} and (|6]) are proportional 
to the molar volume Vm with the same constant of proportionality for all compositions, (a) tests the original shoving model 
(Eq. (HI), according to which the quantity GVm/Tg is the same for all systems, where G is the glass shear modulus. This 
is obeyed to a fairly good approximation, but there is a slight correlation with the glass' Poisson ratio.— This correlation is 
opposite and considerably stronger for the bulk modulus controlled elastic model [(b)] . The energy-landscape analogue of the 
shoving model [Eq. (|6}], which is tested in (c), does not show correlation with the Poisson ratio. 
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